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We present the quantum limits to the magnetic sensitivity of a new kind of magnetometer based 
on biochemical reactions. Radical-ion-pair reactions, the biochemical system underlying the chem- 
ical compass, are shown to offer a new and unique physical realization of a magnetic field sensor 
competitive to modern atomic or condensed matter magnetometers. We elaborate on the quantum 
coherence and entanglement dynamics of this sensor, showing that they provide the physical basis 
for testing our understanding of the fundamental quantum dynamics of radical-ion-pair reactions. 



Quantum physics, in particular quantum information, 
control and measurement are concepts rapidly infiltrating 
biological or biochemical processes. For example, quan- 
tum coherence has been shown to play a fundamental role 
in photosynthesis [H-Q , while quantum measurement dy- 
namics have been demonstrated to underlie radical-ion- 
pair reactions 0-[H| ; the spin-dependent biochemical re- 
actions at the heart of the avian magnetic compass mech- 
anism. Even electron spin entanglement has also been 
addressed with respect to the chemical compass Il2l. 13J , 
paving the way for the dawn of quantum biology [14 1. 

We will here show that spin-selective radical-ion-pair 
reactions are no different than atomic [l5[ or solid state 
quantum sensors [H used in e.g. precision metrology 
[17H19I ]. and in principle are able to offer an exquisite 
magnetic sensitivity. We will establish the fundamental 
quantum limits to the magnetic sensitivity of these bio- 
chemical sensors, and we will elaborate on the fundamen- 
tal decoherence mechanism present. We will show that 
the mechanism damping singlet-triplet coherence also 
damps any electron spin entanglement possibly present, 
as is well understood in precision measurements. In so 
doing, we will also show that the recently appeared en- 
tanglement considerations [12, [l3| are questionable, since 
these works have not taken into account the fundamental 
singlet-triplet decoherence process. Finally, we will com- 
pare the three different and currently competing master 
equations attempting to describe the quantum dynam- 
ics of radical-ion-pair reactions. Based on their ability 
to correctly account for the coherence and entanglement 
dynamics of these reactions, it will be shown that one 
out of the three theories can be eliminated. 



Radical-ion pairs (Fig. 1) [20l - |22j are biomolecular ions 
with two unpaired electrons and any number of magnetic 
nuclei, created by a charge transfer from a photo-excited 
D*A donor-acceptor molecular dyad DA. The magnetic 
nuclei of the donor and acceptor molecules couple to the 
two electrons via the hyperfine interaction, leading to 
singlet-triplet (S-T) mixing, i.e. a coherent oscillation of 
the total electron spin state, also affected by the elec- 
trons' Zeeman interaction with the external magnetic 
field. Singlet-triplet coherence (and its relaxation) is 
studied in a variety of contexts, as e.g. in NMR [23l — 
[25j and quantum dots [26l - [28j . The additional compli- 
cation of radical-ion-pair reaction dynamics is the spin- 
dependent loss of radical-ion pair population due to the 



recombination effect. Charge recombination terminates 
the reaction leading to the neutral reaction products. An- 
gular momentum conservation at this step enforces the 
reaction's spin selectivity: only singlet state radical-ion 
pairs can recombine to reform the neutral DA molecules, 
whereas triplet radical-ion pairs recombine to a differ- 
ent metastable triplet neutral product. If the radical- ion 
pair's lifetime is large enough, the minute Zeeman inter- 
actions can significantly affect the reaction yields, ren- 
dering the reaction a biochemical magnetometer. The 

D*A 



Triplet 

Measurement 
Apparatus 



(DA) T " 




FIG. 1: (Color Online) A donor-acceptor dyad is photoex- 
cited to the state D*A, which after charge separation pro- 



duces the singlet radical-ion pair (D* + A* 



Magnetic in- 



teractions within the molecule induce an S-T interconversion, 
terminated by the recombination event, that results into a 
neutral singlet, DA or triplet product, (DA) T . The funda- 
mental decoherence process is the S-T decoherence, due to the 
coupling to the excited vibrational levels of DA and (DA) T , 
denoted by (DA)* and (DA) T *, respectively. These perform 
two functions: (i) virtual transitions from the radical-ion-pair 
to those levels and back interrupt the unitary magnetic evolu- 
tion, causing S-T decoherence at a rate (ks + kr)/2, and (ii) 
they are a sink of radical-ion-pair population, as real transi- 
tions to them result in the radical-ion pair's singlet (triplet) 
recombination at a rate ks (kr)- 

actual physical signal carrying the magnetic field infor- 
mation can take various forms, for example it can be the 
reaction yield e.g. the singlet, it can be the magnctic- 
field-dependent time evolution of (i) the radical-ion pair 
population, measured e.g in optical absorption experi- 
ments HI, or (ii) the singlet radical-ion-p air p opulation 
measured e.g. in fluorescence experiments [3fJ. No mat- 
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ter what the actual magnetometric signal or the particu- 
lar measurement scheme, the fundamental magnetic sen- 
sitivity of this sensor is determined by general arguments 
no different than similar considerations in atomic sensors, 
i.e. based on the energy resolution during a finite mea- 
surement time, in this case the reaction time. Indeed, 
it is well known [lH [l7l - [l~9j that if we use a single quan- 
tum system to measure an energy splitting E during time 
T, Heisenberg's uncertainty relation (with H = 1) limits 
the precision SE to SE = 1/T. It is also known that 
we can do better than that by simultaneously employ- 
ing N (uncorrelated) quantum systems, leading to an 
improvement by ^J~N . No matter how the measurement 
is performed, the fundamental shot-noise-limited energy 
measurement precision is thus SE = 1/yNT. Now, if 
we let an ensemble of electron spins interact with a mag- 
netic field B for a time T, the limit SE readily trans- 
lates into the limit SB, the smallest measurable change 
of the magnetic field, or equivalently, the magnetic sensi- 
tivity, which is SB = 6E/7, where 7=2.8 MHz/G is the 
electron gyromagnetic ratio. So the fundamental mag- 
netic sensitivity of a system employing N electron spins 
is SB = 1/jVWT. 

Radical-ion pair (RP) reactions start by photo exciting 
the neutral precursors. We assume that at time t = a 
flash of light creates No radical pairs. The RP reaction, 
i.e. the creation of the neutral reaction products, takes 
place at a characteristic time T r , the reaction time. Con- 
sidering for simplicity a mono-exponential decay of the 
radical-ion-pair population of the form N(t) = Noe~ t / Tr , 
the maximum of tyj N(t) (and hence the minimum of 
SB) is realized at t = 2T r , hence, apart from ordcr-of- 
unity factors, the shot-noise-limited magnetic sensitiv- 
ity is SB — l/^\fNoT r . If the photoexcitation-reaction 
cycle is repeated n times for a total measurement time 
t = n(2T r ), averaging the n measurements will yield a 
\fn improvement, thus 



SB = L 
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For an order of magnitude estimate it is noted that con- 
centrations of radical-ion pairs on the order of 10 16 cm -3 
can be obtained, while reaction times up to 10 /is or 
higher are not uncommon, leading, for an active vol- 
ume of 1 cm 3 , to a sensitivity on the order of SB 
10~ 13 / yrfs] G/vTIz. The bandwidth of the magnetome- 
ter will obviously be on the order of l/r. 

All quantum sensors are plagued by decoherence, this 
one being no exception. As explained in detail in [Tlj. 
the fundamental decoherence process of this magnetome- 
ter stems from the continuous measurement of the elec- 
tron spin state performed by the coupling of the radical- 
ion-pair to the reservoir states, which are inherent in the 
molecule DA, as shown in Fig. 1. This coupling results 
in projections on the singlet or triplet subspace, occur- 
ring randomly along the singlet-triplet (S-T) mixing, and 
leading to the damping of S-T coherence. This is an in- 
trinsic, minimum and unavoidable, or in other words fun- 



damental decoherence process, any other extrinsic spin 
relaxation phenomena only adding up and reducing the 
sensitivity (making SB larger). The relevant decoherence 
time turns out to be 0, [TT|] on the order of the reaction 
time T r . This is why the decoherence time does not inde- 
pendently enter (fT]) , in contrast to atomic magnetometers 
[15j . the sensitivity of which is determined by two unre- 
lated time scales, the measurement time and the decoher- 
ence time. Here these two time scales coincide. It is also 
noted that in current biological/biochemical realizations 
the actual sensitivity is far from saturating the limit ((T|) , 
the main reason being the sub-optimal magnetometric 
signal measurement, which is far from being shot-noise- 
limited [2!|[33j]. In that case, i.e. if the measurement's 
signal-to- noise ratio is S/JV < ^/Nq, and setting hence- 
forth t = T r , it follows that 
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The preceding discussion is well-understood in the field 
of precision measurements, for which these quantum lim- 
its are what the impossibility of perpetuum mobile con- 
structs is in the field of thermodynamics. For the skep- 
tic general reader, however, we will give a concrete ex- 
ample elucidating the fundamental impasse posed by 
such limits. One could ask, if the singlet reaction yield 
Ys = Yg(B), clearly dependent on B, is chosen as the 
magnetometric signal, why is the precision SB limited 
by the reaction time? The answer, as mentioned in the 
introduction, is that the shorter the reaction time, the 
less dependent the yield on the magnetic field. This 
can be easily seen, because in the limit T r — > it 
would be Ys — > 1, i.e. the magnetic Hamiltonian has 
no time to mix the electron spin, hence the reaction 
products are all singlet (if the initial state was a sin- 
glet), independently of B. In the general case, and 
assuming equal recombination rates ks = for = k, it 
can be shown [3l| that e.g. the singlet reaction yield 

iB Y S = {l/M)Y^^^\{Qs)nrn\y/{k 2 + O&J. 

where M is the nuclear spin multiplicity, u n are the 
eigenvalues of H with n = 1,2, ...,4M (4 is the multi- 
plicity of the two-electron spin space), w„ m = ui n — u m , 
(Qs)nm is the matrix element of Qs in the eigenbasis of 
% and Qs the 4M-dimensional singlet projection oper- 
ator. Thus, as k — > 00, the sensitivity, dYs/dB, of the 
yield on the magnetic field, encoded in Lu nm , tends to zero 
as 1/k 2 . If the precision in measuring Ys were SYs, then 
the magnetic sensitivity would be SB = SY$/{dYs/dB) 
and would scale as 1/T, 2 , over-satisfying the limit @. 

The limit ^ was derived on the usual assumption of 
uncorrelated single-spin systems. Radical-ion pairs, how- 
ever, are two-electron-spin systems, a possible entangle- 
ment of which could in principle lead to an improvement 
in SB by a factor of \/2 at most, as already known fl7i — 
[i~9j . Neglecting this order-of- unity enhancement, we re- 
peat that ([I]) is a hard limit that cannot be surpassed, 
no matter how ingenious the measurement scheme. Gen- 
eralizing the above example, if we choose to measure an 
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observable, O — O(B), which depends on the magnetic 
field, and the measurement precision of O is SO, then the 
measurement precision of B will be 



SO 



\dO/dB\' 



(3) 



where \dO/dB\ is the sensitivity of the B-dependence of 
the observable O. The achieved limit SBq cannot be 
lower than SB no matter what, i.e. SBq /SB > 1. Re- 
cently, however, in their study of the entanglement dy- 
namics of radical-ion-pair reactions, Briegel and cowork- 
ers [121 ] found a fundamentally different result, severely 
violating the above inequality. These authors introduced 
an observable O — Te, the entanglement lifetime, shown 
to have a steep dependence on the magnetic field B, de- 
picted in Fig. 2b of [12j. We can use this steep de- 



pendence to measure B. To that end we would also need 
the precision ST E of measuring the entanglement lifetime 
T E . Measuring time is like measuring frequency, or en- 
ergy. If we measure a frequency / during T r , the best 
possible precision is 1/T r , again improved by (S/N), i.e. 
Sf = l/T r (S/N). If / 



1/Te, then Sf = 5T E /T E , 



hence ST E = T E /[T r (S/N)}. Using © we find that the 
magnetic sensitivity attained through a measurement of 
T E is 



5B Te = 



1 



T 2 



1 



(S/N) T r \dT E /dB\' 



(4) 



It must be SBt e > SB, thus the inequality \dT E /dB\ < 
jT E must be satisfied, otherwise the fundamental limit 
([2]) is violated. However, in the calculation of the entan- 
glement lifetime's dependence on B found in [l2j], there 
is nothing limiting the slope \dT E /dB\, i.e. by zoom- 
ing into the region of the steep B-dependence, the slope 
increases arbitrarily, and hence the violation of ([2]) is 
arbitrarily high. We will now clarify the root of this 
inconsistency. Briegel and coworkers considered T E to 
be determined just by the magnetic interactions present 
in the radical-ion-pair. This is an erroneous association 
that disregards S-T decoherence inherent in radical-ion 
pairs. As in all precision measurements, the spin coher- 
ence time is the upper limit of the entanglement life- 
time fl5l - F]~9| . Parenthetically, this is the reason why en- 
tangled states are hardly helpful in surpassing the stan- 
dard quantum limits in quantum parameter estimation 
[HI HU, i.e. any entanglement decays at least as fast as 
the coherent superpositions. In contrast, the authors in 
fl2j calculate an entanglement lifetime that bears no re- 
lation to the relevant decoherence time, which is on the 
order of T r , in fact they find that in cases T E ^> T r . To 
further elucidate the root of the problem, we will con- 
sider a very simple example of a radical-ion pair with- 
out nuclear spins, just the two electrons, starting out in 
the singlet state |^o) = \S) = (| + -) - | - +})/v / 2- 
We suppose that the only magnetic interactions are the 
Zeeman interactions of the two electrons, and in order 
to induce S-T mixing, we consider the Ag mechanism 
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FIG. 2: (Color Online) (a) Time evolution of (Qs)t and £t 
for the case of no reaction, ks — kr = 0, calculated from 
dp/dt = —i[H,p]. (b) Time evolution of (Qs)t and (c) St 
for the realistic case of non-zero recombination rates ks = 
kr = Slj, calculated with the full master equation dp/dt = 
—i[H, p] + C(p) of the Jones-Hore theory [Io| and the Kominis 
theory [ill ] . The full master equation of the traditional theory 
[37], HH leads to the same result as in case (a). In the right y- 
axis of (c) we plot the normalization of p, in order to elucidate 
the fact that there still exists a substantial number of radial- 
ion pairs at the time when the predictions of the two theories 
(Jones-Hore and Kominis) start to significantly deviate from 
the third (traditional). 



|32j, i.e. the Larmor frequencies of the two electrons are 
taken to be slightly different. Thus the magnetic Hamil- 
tonian is simply H = uj\Si z + W2«2z- If we completely 
disregard the reaction dynamics and for the moment as- 
sume that the evolution of the density matrix is driven 
solely by T-L, we find the the initial state |^q) evolves to 
\i/j t ) = e -^*/2[(i + e^ 4 )^) + (1 - e iSut )\T )]/2, where 
\Tq) = (| H — ) + I — h))/\/2 is the zero-projection state of 
the triplet manifold spanned by \T ) and \T±) = \ ± ±). 
As expected, H induces a coherent oscillation S <^ T at 
a frequency Suj = u>i — uj^. It is easily shown that the 
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concurrence [Hj], i.e. the overlap of \ipt) with the time- 
inverted \tp t ) — cr y|V't) is at all times C t — \(ipt\i>t) \ = 1- 
Thus the electron spin state starts out and remains max- 
imally entangled and maximally coherent. In reality, 
however, radical-ion pairs suffer a continuous loss of S-T 
coherence due to the continuous quantum measurement 
induced by the recombination dynamics. In a single- 
molecule picture, this intramolecule measurement results 
in random quantum jumps either to the singlet or the 
triplet manifold of the RP spin states. In an ensemble of 
RPs this is equivalent to an improper mixture of singlet 
and triplet states, with the concomitant loss of S-T co- 
herence and entanglement. The above considerations are 
illustrated in Fig. [SJ 

In Fig. [2^, we plot the expectation value of Qs and the 
entanglement £ t — £(C t (pj) [36| for a fictitious radical- 
ion pair, for which k$ — &t = 0, i.e. when the evolution 
of p is driven just by the magnetic Hamiltonian % and the 
singlet and triplet recombination channels simply do not 
exist. As explained above, the expectation value (Qs)t 
exhibits an undamped S-T oscillation with frequency Su, 
while the entanglement remains at its maximum value of 
one indefinitely, i.e. both the coherence and the entan- 
glement lifetimes seem to be infinite. In Figs. and 
|2J; we consider the physical case, i.e. we turn on the re- 
combination channels and take k$ — &t = Suj. Before 
explaining the result, we note that there currently exist 
three different theories attempting to describe the funda- 
mental dynamics of radical-ion-pair reactions, the tradi- 
tional theory of spin chemistry derived in (37l . l38j , a new 
theory based on quantum measurement considerations 
S El and a theory based on a different interpretation of 
the quantum measurement introduced in [lOj. All three 
theories are of the form dp/dt = —i[H, p] + C(p), where p 
is the electron and nuclear spin density matrix and C(p) 
the super-operator describing the effect of the reaction. 
The three theories differ in the particular form of C(p), 
while they obviously agree in the imaginary scenario of 
no reaction, since when k$ — kx — it is C(p) = 0, 
and their common result was depicted in Fig. [2^,. Turn- 
ing on the reaction channels, however, vast differences 
show up. In Figs. [2Jd and it is seen that both S-T 
coherence and entanglement, respectively, decay in both 
the Jones-Hore and the Kominis theory, (albeit on a dif- 



ferent time scale for each theory), exactly because the 
intra-molecule measurement dynamics transform coher- 
ent superpositions into incoherent mixtures. In contrast, 
the traditional theory leads to exactly the same result 
as in Fig. [2^, the case of no reaction! The reason is 
that the traditional theory is founded on the erroneous 
physical assumption that nothing happens to surviving 
(unrecombined) radical-ion pairs. Indeed, since in the 
considered case of fcs,fcr ^ all radical- ion pairs even- 
tually recombine away, i.e. Tr{p} — > as t — > oo, we 
normalize p by Tr{p}, in order to describe the entan- 
glement and singlet expectation value, as a function of 
time t, of those molecules that have not recombined un- 
til time t. What is seen in Fig. [5]d is that the oscillations 
°f (Qs)t are damped away while (Qs}t decays towards 
(Os)oo, where (Q s )t = Tr{Q s p}/Tr{p}. Similarly (Fig. 

£ t = £(C*(p/Tr{/9})) also decays to zero. Thus in re- 
ality, both the coherence and the entanglement lifetimes 
are finite and on the order of the reaction time. 

Put differently, the premise of the entanglement life- 
time calculations of [12], a premise embedded in the tra- 
ditional theory, is that until they recombine, RPs evolve 
unitarily under the action of W only. In contrast, the 
other two theories explicitly address the fact that the re- 
combination dynamics are not only responsible for trans- 
forming RPs into neutral products, but also for inducing 
S-T decoherence of RPs until they recombine. The ab- 
sence of this physical mechanism from the fundamental 
theoretical description of RP reactions leads to unphys- 
ical coherence and entanglement dynamics, and in their 
turn, these lead to violation of fundamental quantum lim- 
its. Furthermore, the entanglement considerations of Ve- 
dral and coworkers [r| are also based on the traditional 
theory, as the authors themselves state, rendering their 
qualitative and quantitative conclusions questionable. 

In summary, we have set the entanglement dynamics 
of the chemical compass on a fundamentally sound ba- 
sis. In so doing, we have shown that one out of the 
three currently competing master equations purporting 
to describe the quantum dynamics of radical-ion-pair re- 
actions cannot stand as a fundamental theory, since it 
violates well-established fundamental quantum limits of 
precision measurements. 
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